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SUPPRESSION OF BUOYANCY IN GASEOUS MEDIA 
AT HIGH TEMPERATURES 


Suleyman A. Gokoglu and Maria A. Kuczmarski 
National Aeronautics and Space Administration 
Glenn Research Center 
Cleveland, Ohio 44135 

Introduction 

Natural convection, induced by buoyancy forces in the presence of a gravitational field, is 
driven by density gradients generated by either temperature or concentration gradients. If the 
gradients are perpendicular to the gravity vector, the fluid is unconditionally unstable. Above a 
horizontal surface, however, the transport of heat and mass in a quiescent fluid will be governed 
solely by conduction and diffusion until buoyancy forces become strong enough to initiate 
motion. In our discussions below, we limit our considerations only to thermal gradient cases 
where the heated surface is below the colder fluid and is horizontal with respect to the gravity 
vector. The onset and development of buoyant flow under such conditions depends quantitatively 
on Rayleigh number, defined as 


Po M 


(1) 


where p = p 0 - p. If the Rayleigh number is smaller than a critical Rayleigh number Rac, then 
the fluid remains still (known as the Rayleigh-Benard criterion). 

The Rayleigh number can also be expressed as a product of the Grashof number and the 
Prandtl number, respectively defined as 


Gr = 


A p gL 3 

P 0 v 2 



P c P 

X 


( 2 ) 


such that Ra = Gr • Pr. Under small normalized temperature differences at constant pressure, 
Boussinesq conditions apply (Boussinesq, 1903); i.e., it is commonly understood that the 
governing conservation equations consist of the following approximations: (1) density is 
assumed constant except when it directly causes buoyant forces, (2) all other fluid properties are 
assumed constant, and (3) viscous dissipation is assumed negligible, as was first used by 
Oberbeck (1879). The normalized density difference can then be approximated by the first term 
of its Taylor-series expansion; i.e., 

— * P t (T-T q )= J3 t AT (3) 

Po 

where (3j, defined as (3 t = - (l/p)*(3p/3T) p , is the coefficient of thermal expansion at constant 
pressure. If a temperature gradient (i.e., density stratification) is too large for diffusion to sustain, 
such that a Rac is exceeded, then finger-like intrusions of the lighter gas into the heavier gas, and 
the heavier gas into the lighter gas, result rapidly. 
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The validity of the Boussinesq approximation (Gray and Giorgini, 1976) and departures from 
the Boussinesq approximation (Paolucci and Chenoweth, 1987) have been studied in the 
literature. These studies led to the specification of explicit conditions under which the neglect of 
various terms in the traditional Boussinesq approximation can be justified. 

Goldstein and Volino (1995) give an extensive description of various investigations in the 
literature for the onset of convection and flow development in horizontal layers. Many questions 
are raised due to discrepancies among different flow structures. Some studies had free surfaces, 
while others used rigid surfaces. The heated surface was applied either a step change or a linear 
increase in temperature, or a step change in heat flux. Some of the investigators report plunging 
sheets, while others observe isolated plumes and thermals, and still others see both. They use 
different criteria and different experimental observations for the basis of determining the onset of 
instability. Several authors report widely ranging values of Ra c , easily from 1000 to 3000. 


Temperature-Dependent Properties for Gases and Its Consequences on Buoyancy 

Albeit the uncertainty on Rac, even under simpler Boussinesq conditions, the problem of 
determining a critical Rayleigh number for non-Boussinesq conditions ( large normalized 
temperature differences), using high-order eigenvalue formulations/solutions, is quite intriguing 
and remains to be thoroughly treated. Upon extension to large temperature differences, the 
temperature dependence of the thermodynamic and transport properties becomes a matter of 
significance and can no longer be neglected. 

More recently, Frohlich et al. (1992) investigated large departures from the Boussinesq 
approximation in the Rayleigh-Benard problem. They determined that temperature-dependent 
properties affect the mid-plane symmetry and lead to quantitative and qualitative changes of the 
flow, in particular, in the bifurcation from the conduction state. Their most important finding 
states that the non-Boussinesq behavior becomes more confined to small regions in the flow, and 
hence, the relative importance of these effects is expected to diminish for higher Rayleigh 
numbers. However, their study fails to determine that gaseous flows tend to become more stable 
towards the conduction state at larger normalized temperature differences, due primarily to the 
dampening effect of gas viscosity at higher temperatures, as investigated in this paper. 

Here, we limit our considerations only to gaseous media because the temperature dependence 
of thermodynamic and transport properties for gaseous systems leads to some interesting 
consequences in regards to flow stability and onset of buoyant convection. For ideal gases, it is 
readily shown that (3 t = 1/T. The power-law approximations of the temperature dependence of 
other gaseous properties are typically taken as (Svehla, 1962): p ~ T 1 , p ~ T 0 ' 7 , X ~ T 0 ' 85 , c p ~ 
I 0 ' 19 , D ~ T 1 ' 7 . These powers are also consistent with the kinetic theory of gases (Kennard, 
1938; Hirschfelder, Curtiss, and Bird, 1954), and are accurate descriptions for the properties of 
air in the 600 to 1600 K range (Brown and Donoughe, 1951; Kays and Crawford, 1980). 

Interestingly, the dimensionless numbers based on the relative importance of momentum, 
heat, and species mass diffusion, namely the Prandtl (Pr), Schmidt (Sc=v/D=p/pD), and Lewis 
numbers (Le=ot/D=X/pcpD), become temperature independent, i.e., the power of temperature is 
smaller than 0.04. Hence, ignoring variable property effects for many convection problems is 
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justified. Since the temperature dependences of the Rayleigh and Grashof numbers become 
about the same, we focus our attention to Gr from here on. 


For large normalized temperature differences, the variable (temperature-dependent) gas 
properties are typically accounted for by evaluating their values at a reference temperature T*, 
which is customarily taken as the arithmetic mean of the hot and cold surface temperatures in the 
system (Kays and Crawford, 1980; Bird, Stewart, and Lightfoot, 1963); i.e., 


t . = r, +0.5(7; -rj 


2 


(4) 


Paolucci and Chenoweth (1987) also found in their calculations that this reference temperature 
describes average properties more accurately than other alternatives. The temperature 
dependence of the Grashof number can now be expressed as 


Gr = /3 t AT 



AT gL' 
T. v 2 


AT 1 

rp np 3.4 

J- * J- * 


= AT 


^—4.4 


(5) 


It is indeed noteworthy that the temperature dependence of v, the gas kinematic viscosity, is 
primarily responsible for making the Grashof number so strongly temperature dependent. 

AT 

Typically, the Boussinesq approximation is valid if <0.1 (6) 


For cases where T 0 is room temperature around 300 K, AT is thus limited to smaller than 35 K. 
For most heat transfer applications, AT is much larger than 35 K, and for combustion, it will 
routinely exceed 1000 K. Hence, one has to be cautious in the use of such things as the Rayleigh 
stability criteria or heat transfer correlations based on the Rayleigh (Grashof) number. 

For quantifying the temperature dependence of the Grashof number, consider a constant- 
property Grashof number, Gr„, for which the gas properties are evaluated at the cold-surface 
temperature, T 0> and the normalized temperature difference is taken as unity (or AT = 1 K); i.e.. 


Gr 0 = 


j _gV_ 

t v: 


(7) 


Then the variation of the Grashof number with temperature at larger normalized temperature 
differences can be quantified by normalizing the Grashof number with respect to Gr 0 . Thus, the 
normalized Grashof number, defined as GR, has the following temperature dependence: 


Gr 

GR = = AT 

Gr o 



V- 1 * J 


( 8 ) 


Note that Eq. 8 is independent of the type of gas and the pressure of the system as long as the 
power-law descriptions for the temperature dependence of gas properties are adequate. Figure 1 
below describes the behavior of GR for three different cold surface temperatures. The main 
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curve is the one labeled for cold-surface temperature of T 0 = 300 K and forms the basis of most of 
the discussions presented in this paper. Other curves labeled at lower cold-surface temperatures 
apply to arguments related to cryogenic applications. 

Figure 1 has a number of interesting implications. It suggests that if a system becomes 
unstable when the hot-surface temperature is increased, then there are two distinct critical hot- 
surface temperatures for the onset of buoyant convection, one at a lower temperature and another 
at a higher temperature. Below and above these temperatures, the system is stable, i.e., 
conduction-dominated. It seems that, for a fixed cold-wall temperature, the lower the critical 
low-temperature is for a system, the higher the critical high-temperature is. The question of 
whether there is any significant difference in the critical Rayleigh numbers for the onset of 
convection at the lower and higher temperature regions will be explored below. 
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Figure 1. — Variation of normalized Grashof number with hot surface temperature 
(eq. 8) describing the effect of temperature-dependent gas properties. 

Numerical Modeling Results 

In order to provide an initial check for the validity of fig. 1, we have built a 2-D model using 
the computational fluid dynamics code FLUENT. The rectangular geometry and the conditions 
are taken to represent the traditional Benard problem, where the distance between the bottom hot 
and the top cold surfaces is fixed at 1.25 cm and the aspect ratio is chosen as 6. The vertical 
side-walls are insulated (adiabatic), the top cold-surface temperature is fixed at 300 K, and the 
bottom hot-surface temperature is varied parametrically at fixed levels. We use either neon or 
argon whose temperature-dependent properties are fully accounted for by using curve fits. 
Pressure is also fixed, but parametrically varied to explore, under steady-state conditions, 
whether the two different regimes are dominated either by conduction or buoyant convection. 

Our numerical experiments for Ne under atmospheric conditions reveal that Ne is always 
stable for all hot-surface temperatures. This finding is not surprising since the largest Rayleigh 
number for Ne at the peak temperature of -476 K of figure 1 is calculated to be -1500. This is 
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smaller than Ra c values reported by Frohlich et al. (1992) and Goldstein and Volino (1995), 
which are larger than 1800 and 2000, respectively. These results demonstrate that one cannot 
make the Ne system unstable to onset buoyant convection by simply raising the hot-surface 
temperature. However, when we switch the working gas to Ar, the system is unstable at the hot- 
surface temperature of 476 K, as expected and depicted in fig. 2. In all of our calculations where 
the gas becomes unstable, recirculating cells with aspect ratios of about unity are generated, 
creating warm and cool fingers. The stable, quiescent cases, however, result in almost linear 
temperature profiles, where the slight departures from linearity are due to temperature-dependent 
gas properties. Note that the color scale is linear and applies the same way to all figures. 


(a) Ne temperature field; 
stable, quiescent gas. 



Figure 2. — Ne and Ar temperature and velocity fields for T h = 476 K (peak temperature for fig. 1) 
and p = 1 atm. The color scale is linear. 

Based on the molecular-weight ratio of about 2 between Ar and Ne, the dependence of the 
Grashof number on the square of density and that Ne has smaller viscosity than Ar, we expect Ne 
to be unstable for the same hot- surface temperature of 476 K, but at 5 atm. Indeed, figure 3 
demonstrates the resulting slightly stronger velocity field with its associated temperature field. 



(a) Ne velocity field, 
V m ax= 12.0 cm/s. 


(b) Ne temperature field. 


Figure 3. — Ne velocity and temperature fields for T h = 476 K and p = 5 atm. 

In the next set of calculations, we successively increase the pressure for Ne gas 
from atmospheric pressure, for which Ne stays stable for all hot-surface temperatures, to as high 
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as 5 atm. The results are schematically shown in fig. 4 below and reproduce fig. 1 in reversed 
shape. Indeed, the implications of fig. 1 are numerically demonstrated again. Furthermore, one 
can readily observe from fig. 4 that some systems remain stable with respect to the onset of 
convection regardless of how large the normalized temperature difference becomes. In fact, 
beyond a certain maximum hot-surface temperature, the tendency of systems to become more 
stable towards a conduction-dominated state increases as the hot-surface temperature increases. 



Figure 4. — A schematic of calcu- 
lated stable and unstable regions, 
and hot-surface temperatures at 
which transition occurs for Ne at 
various pressure levels. 


The proper determination of the critical Rayleigh numbers for the onset of convection at the 
lower and higher temperature regions requires transient analysis. Such calculations (subject of 
our current work) can also quantify the required heating rates of the hot-surface in order to pass 
through the unstable regime in the middle to reach the stable regime at higher temperatures. We 
tabulated below the calculated Rayleigh numbers corresponding to the transition temperatures 
depicted in fig. 4, cognizant of the limitations of using steady-state analysis. 


P, atm 

T h ,K 

Ra c 

1.2 

372 

761 

1683 

1706 

1.5 

338 

1010 

1812 

1672 

3.0 

308 

1892 

1853 

1613 

5.0 

303 

2791 

1996 

1570 


It is noteworthy that, for the same pressure pairs, the critical Rayleigh numbers depart more 
from one another as the pressure increases; i.e., as the T h differences increase. This may be due 
to the fact that there are larger inaccuracies in the calculation of Rac at larger differences between 
hot- and cold-wall temperatures, since properties evaluated at the average temperature are used. 
In fact, the Boussinesq approximation is valid in only two of these cases (see eq. 6). It should 
also be noted, however, that all critical Rayleigh values are in the range of those most recently 
reported by Goldstein and Volino (1995), which vary depending on the measurement technique, 
and are given as 2673 on the basis of thermocouple, 2700 on the basis of liquid crystal, and 1265 
on the basis of interferometry measurements. 
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Naturally, the level of convection is weaker at the hot-surface temperature where transition 
from stable to unstable occurs (determined numerically to within 1 K). An example of such 
small velocities is shown below in fig. 5 relative to those shown in figs. 2 and 3. We determined 
this transition temperature of 308 K for Ne at 3 atm, such that at Th= 307 K the temperature 
profile was completely conduction dominated. 
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(a) Ne velocity field, 
Vmax = 0.5 cm/s. 



(b) Ne temperature field. 


Figure 5. — Ne velocity and temperature fields for the transition Th= 308 K and p = 3 atm. 


Some Practical Implications 

Generalization of the Rayleigh-Benard criterion for the onset of convective instability to non- 
Boussinesq conditions is long overdue in the heat transfer/fluid physics community. We list 
below a few examples where such information is not only relevant, but critical. 

The application of existing heat transfer correlations involving the Rayleigh and Grashof 
numbers for engineering purposes is a routine practice, even under conditions where the 
Boussinesq approximation is no longer valid. Though for liquid systems there is little concern, 
we may not be so fortunate for gaseous systems, especially if the hot-surface temperature exceeds 
-500 K. According to fig. 1, we predict that the tendency for buoyant convection will be 
suppressed more at higher temperatures. If a practicing engineer uses recommended traditional 
methods to calculate the relevant Rayleigh (Grashof) numbers for hot-surface temperatures much 
beyond 500 K, he/she may be misled to a wrong conclusion, because these numbers will be larger 
at even lower temperatures through which the system will have to pass to reach the higher levels. 
One may conclude a stable, quiescent, conduction-dominated system behavior based on the 
calculated smaller numbers, whereas the system may become unstable and convection-governed 
while heating up to the higher temperature. Similarly, even if both predictions and reality may 
result in the unstable, convection-governed regime, the intensity of buoyant convection may be 
totally misjudged by not considering the history of the system during heat-up. The consequences 
of such misjudgments can be severe in both earth and space applications with partial levels of 
gravitational acceleration, such as Lunar or Martian. 

For materials processing involving gaseous media, particularly chemical vapor deposition 
(CVD), buoyancy-driven convection is avoided to improve deposition rate uniformity and 
deposit quality. The effect of temperature-dependent gas properties on flow dynamics in the 
CVD process has long been recognized (Rosenberger, 1987). In horizontal CVD reactors, 
buoyancy leads to a helical flow structure. In vertical, impinging-jet or stagnation-flow CVD 
reactors, buoyancy is suppressed by either enhanced reactant flow rates or by lowering the 
pressure, both of which result in reduced efficiencies (Gokoglu, 1992). General perception in 
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CVD is that if the temperature of the hot-substrate increases, the tendency for buoyant 
convection, which is undesirable, also increases. Yet, our analysis above suggests that if the hot- 
substrate is heated to even higher temperatures (in excess of about 500 K), this may even be 
beneficial in suppressing buoyancy. The challenge would be to heat the substrate to these high 
temperatures rapidly enough to avoid the development of buoyant convective instabilities at 
lower temperatures during the ramp-up. This should be possible with currently available 
technology used in rapid thermal processing (RTP) of materials. Ironically, because high 
substrate temperatures are sometimes harmful to the desired deposit quality, especially for 
electronic CVD materials, the low-temperature CVD processes employed to circumvent this 
problem may be more prone to the buoyant convection problem. For applications of CVD to the 
growth or coating of structural materials, where the use of high substrate temperatures is not as 
critical an issue, the results presented here are directly relevant. 

There are serious consequences of our proposed research for partial-gravity combustion 
applications in space and for the combustion community at large. Diffusive-transport-controlled 
flames with planar symmetry create environments with large normalized temperature differences 
where the traditional application of Rayleigh-Benard criterion for the onset of convective 
instability is quite inadequate (Baumstein and Fendell, 1998; Fendell and Mitchell, 1996). Upon 
extension to non-Boussinesq conditions, it may prove necessary to regard the critical Rayleigh 
number as a function of the Prandtl and Lewis numbers, as well as the normalized temperature 
difference between the diffusion flame and the cold upper surface. The thermodynamic 
dependence of the transport properties becomes a matter of import. 

Applying traditional scaling approaches by the use of dimensionless numbers to micro-scale 
combustion has been shown to be difficult (Sitzki et al., 2001). Critical dimensions of such 
micro-combustors and the heat loss mechanisms associated for optimal operation require a more 
careful analysis. The results presented here suggest that designers of micro-combustors select the 
temperatures of operation and dimensions between hot and cold surfaces more judiciously in 
order to better control conduction versus convection-dominated heat transfer. 

Figure 1 also suggests that the behavior of gaseous systems under cryogenic conditions is 
quite different for the onset of instability. For these systems, where room temperature may be the 
“hot” surface temperature and cryogenically cooled sensors or electronic components may 
constitute the cold surfaces, heat transfer mechanisms may be quite different than those expected 
from simpler, Boussinesq type analyses, because variable gas properties seem to dominate the 
flow behavior. Indeed, for the curve labeled as T 0 = 50 K, the system is drastically dampened by 
the higher viscosity of the gas when the hot surface temperature Th is only 300 K! If the desired 
mode of heat transfer is conductive, then the designers of such devices may highly benefit from 
the results presented here by utilizing even larger temperature differences. 


Concluding Remarks 

Our 2-D numerical work utilizing the computational fluid dynamics code FLUENT 
demonstrate that modeling can provide insightful information to understand Rayleigh instability 
at large temperature differences in gaseous media. The steady state analysis is believed to be 
sufficient for describing the phenomena, though not as accurate as transient calculations, which 
we will report on soon. We show here that under normalized temperature differences larger than 
justifiable for the Boussinesq approximation, there are terrestrial situations where a system 
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cannot ever be made unstable with respect to the onset on buoyant convection no matter how 
large the temperature (density) difference becomes at a given pressure! This behavior is due to 
temperature-dependent gas properties. Highly temperature-sensitive kinematic viscosity is 
primarily responsible for counteracting the tendency toward instability and dampens convection 
by making the gas more “viscous” at higher temperatures. This compensation of the buoyant 
force by the viscous force exhibits itself by the formation of a peak hot-surface temperature 
beyond which a system will tend to be more stable as the hot-surface temperature increases. We 
also show, for a given pressure, that there are two distinct critical hot-surface temperatures for the 
onset of buoyant convection, one at a lower temperature and another at a higher temperature, 
below and above which, respectively, the system is stable, i.e., conduction-dominated. 

Our current work is also exploring the behavior in sealed enclosures, where the pressure is 
allowed to rise during the transient heating of the hot surface. The rate of pressure rise depends 
on the heat-loss rate of the system relative to the heat-up rate of the hot surface and the diffusion 
rate in the gas mixture. While the higher temperatures help dampen the buoyancy effect, the 
higher pressure oppose it. The effect of aspect ratio naturally becomes another parameter in this 
regard. 


Nomenclature 

c p constant-pressure heat capacity 

D species mass diffusivity 

g acceleration of gravity 

Gr Grashof number, Eq. (2) 

L characteristic length (height) between hot and cold surfaces 

Le Lewis number (cx/D = thermal diffusivity / species mass diffusivity) 

M molecular weight 

p pressure 

Pr Prandtl number, (v/oc = momentum diffusivity / thermal diffusivity) 

Ra Rayleigh number, Eq. (1) 

Sc Schmidt number (v/D = momentum diffusivity / species mass diffusivity) 
T temperature 

V velocity 

a thermal diffusivity, X/(pc p ) 

(3t coefficient of thermal expansion at constant pressure 
X thermal conductivity 

p dynamic viscosity 

v kinematic viscosity (v = p/p) 

p density 

Subscripts 

c cold surface condition 

h heated, hot surface condition 

o initial ambient or cold surface condition 

p constant pressure condition 

* average, arithmetic mean temperature condition 


NAS A/TM— 2003-2 1 2695 


9 



References 


Baumstein, A. and Fendell, F. (1998). Strain-Rate-Free Diffusion Flames: Initiation, 
Properties, and Quenching. Combust. Sci. and Tech., Vol. 132, 157-198. 

Bird, R.B., Stewart, W., and Lightfoot, E.N. (1963). Transport Phenomena. John Wiley, 
New York. 

Boussinesq, J. (1903). Theorie Analytitique de la Chaleur, Vol. 2, Gauthier-Villars, Paris. 

Brown, W.B. and Donoughe, P.L. (1951). Tables of Exact Laminar Boundary Layer 
Solutions When the Wall is Porous and Fluid Properties are Variable. NACA TN 2479. 

Fendell, F.E. and Mitchell, J.A. (1996). Feasibility of Planar, Unidirectional-Flow Diffusion 
Flames in Earth Gravity. Combust. Sci. and Tech. Vol. 120, 83-117. 

Frohlich, J., Laure, P., and Peyret, R. (1992) Large Departures from Boussinesq 
Approximation in the Rayleigh-Benard Problem. Phys. Fluids A 4 (7), 1355-1372. 

Gokoglu, S.A. (1992). Chemical Vapor Deposition Modeling for High Temperature 
Materials. Chemical Vapor Deposition of Refractory Metals and Ceramics II, Eds. Besmann, 
T.M., Gallois, B.M. and Warren, J.W., Mater. Res. Soc. Symp. Proc. 250 , pp. 17-28. 

Goldstein, R.J. and Volino, RJ. (1995). Onset and Development of Natural Convection 
above a Suddenly Heated Horizontal Surface. Journal of Heat Transfer 117 , 808-821. 

Gray, D.D. and Giorgini, A. (1976). The Validity of the Boussinesq Approximation for 
Liquids and Gases. Int, J. Heat Mass Transfer, Vol. 19, 545-551. 

Hirschfelder, J.P., Curtiss, C.F., and Bird, R.B. (1954). Molecular Theory of Gases and 
Liquids, J. Wiley and Sons, New York, NY. 

Kays, W.M. and Crawford, M.E. (1980). Convective Heat and Mass Transfer. McGraw Hill, 
2 nd Edition, New York. 

Kennard, E.M. (1938). Kinetic Theory of Gases. McGraw Hill, New York. 

Ng, D. (1999). Temperature Measurement of a Miniature Ceramic Heater in the Presence of 
an Extended Interfering Background Radiation Source Using a Multiwavelength Pyrometer. 
NASA/TM — 1999-208854. 

Oberbeck, A. (1879). Uber die Warmeleitung der Flussigkeiten bei Beriicksichtigung der 
Stromungen infloge von Temperatur Differenzen. Ann. Phys. Chem. 7 , 271-292. 

Paolucci, S. and Chenoweth, D.R. (1987). Departures from the Boussinesq Approximation in 
Laminar Benard Convection. Phys. Fluids 30 (5), 1561-1564. 

Rosenberger, F. (1987). Flow Dynamics and Modelling of CVD. Proceedings of the Tenth 
Int. Conf. On Chemical Vapor Deposition, Ed. G. W. Cullen, The Electrochemical Soc., 
Pennington, NJ, 1 1-22. 

Sitzki, L, Borer, K., Schuster, E., Ronney, P.D., Wussow, S. (2001). Combustion in 
Microscale Heat-Recirculating Burners. Third Asia-Pacific Conference on Combustion, Seoul, 
Korea. 

Svehla, R.A. (1962). Estimated Viscosities and Thermal Conductivities of Gases at High 
Temperatures. NASA Technical Report R-132. 

Weinbaum, S. (1964). Natural Convection in a Horizontal Circular Cylinder. J. Fluid Mech., 
Vol. 18, 409-137. 


NAS A/TM— 2003-2 1 2695 


10 



REPORT DOCUMENTATION PAGE 

Form Approved 
OMB No. 0704-0188 

Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, 
gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson 
Davis Highway, Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-0188), Washington, DC 20503. 

1. AGENCY USE ONLY ( Leave blank) 

2. REPORT DATE 

November 2003 

3. REPORT TYPE AND DATES COVERED 

Technical Memorandum 

4. TITLE AND SUBTITLE 

5. FUNDING NUMBERS 


Suppression of Buoyancy in Gaseous Media at High Temperatures 


6. AUTHOR(S) 

Suleyman A. Gokoglu and Maria A. Kuczmarski 


WB S-22- 1 0 1-42-09 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
John H. Glenn Research Center at Lewis Field 
Cleveland, Ohio 44135-3191 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 


E-14197 


9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Washington, DC 20546-0001 


10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 


NASA TM— 2003-212695 


11. SUPPLEMENTARY NOTES 


Prepared for the Fourth International Symposium on Scale Modeling (ISSM-IV) sponsored by the National Center 
for Microgravity Research on Fluids and Combustion, Cleveland, Ohio, September 17-19, 2003. Responsible person, 
Suleyman A. Gokoglu, organization code 6711, 216-433-5499. 


12a. DISTRIBUTION/AVAILABILITY STATEMENT 

Unclassified - Unlimited 

Subject Category: 34 Distribution: Nonstandard 

Available electronically at http://gltrs.grc.nasa.gov 

This publication is available from the NASA Center for AeroSpace Information, 301-621-0390. 


12b. DISTRIBUTION CODE 


1 3. ABSTRACT (Maximum 200 words) 

Consider a rectangular box filled with a fluid having a heated bottom and a cold top surface, and insulated side- walls 
(Benard problem). As the temperature difference between the horizontal top and bottom surfaces increases, a critical 
condition, defined quantitatively by the Rayleigh number, is reached beyond which density stratification can no longer 
be sustained by conduction and the fluid disrupts from its stable, quiescent state into an unstable, convective mode in 
which lighter and heavier gas mix. This paper suggests that such a statement is not necessarily true for gaseous media 
under normalized temperature differences that are much larger than justifiable for the Boussinesq approximation! In fact, 
there may be situations where a system cannot ever be made unstable with respect to the onset on buoyant convection no 
matter how large the temperature (density) difference becomes at a given pressure even under normal gravity ! This 
“unexpected” behavior is primarily attributed to highly temperature- sensitive kinematic viscosity which counteracts the 
tendency toward instability and dampens convection by making the gas more “viscous” at higher temperatures. This 
compensation of the buoyant force by the viscous force exhibits itself by the formation of a peak hot-surface temperature 
beyond which a system will tend to be more stable as the hot- surface temperature increases. 


14. SUBJECT TERMS 

Taylor instability; Stability; Buoyancy-driven flow; Buoyancy 


15. NUMBER OF PAGES 

16 


16. PRICE CODE 


17. SECURITY CLASSIFICATION 
OF REPORT 

Unclassified 


18. SECURITY CLASSIFICATION 
OF THIS PAGE 

Unclassified 


19. SECURITY CLASSIFICATION 
OF ABSTRACT 

Unclassified 


20. LIMITATION OF ABSTRACT 


NSN 7540-01-280-5500 


Standard Form 298 (Rev. 2-89) 
Prescribed by ANSI Std. Z39-18 
298-102 


